Ghost excitonic insulator transition in layered graphite 
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Some unusual properties of layered graphite, including a linear energy dependence of the quasi- 
particle damping and weak ferromagnetism at low doping, are explained as a result of the proximity 
of a single graphene sheet to the excitonic insulator phase which can be further stabilized in a doped 
system of many layers stacked in the staggered (ABAB . . .) configuration. 



o 
o 

(N 
o 

Q 



> 
o 



C3 



o 
o 



The continuing interest in carbon-based materials has 
intensified both experimental and theoretical efforts to 
understand their electronic properties and interaction- 
driven transitions in these systems. 

The problem of electronic instabilities in a single sheet 
of graphite has been previously studied by means of the 
Hartree-Fock and density functional methods with a fo- 
cus on the short-ranged (including on-site and nearest 
neighbor) Hubbard- like repulsive interactions, and the 
inter-site repulsion was found to favor a charge density 
wave (CDW) ground state . This analysis did not pro- 
vide, however, a proper account of the Coulomb forces 
that remain long-ranged due to the lack of conventional 
screening in semimetals, such as graphite. 

In contrast, the authors of Ref. specifically focused 
on the role of the Coulomb interactions. From the renor- 
malization group calculation in the second order in the 
dimensionless Coulomb coupling g = 2ite 2 /eqv they con- 
cluded that the renormalized coupling monotonously de- 
creases at low energies and, therefore, it can not cause 
any instability of the gapless paramagnetic ground state 
of graphite. 

Besides, the authors of Ref. Q suggested a possi- 
ble explanation for the experimentally observed || lin- 
ear energy dependence of the quasiparticle damping (de- 
fined as the imaginary part of the quasiparticle disper- 
sion e — E p + iT,(e) ) which they found to behave as 
E(e) ~ e/ln 2 e at low energies. However, the large esti- 
mated value of the bare coupling constant (g > 10) calls 
these results into question and warrants further investi- 
gation. 

In the present Letter, we revisit the problem of the 
Coulomb interacting electrons in layered graphite and 
study the nature of the ground state and the quasipar- 
ticle spectrum at strong coupling. This time around, we 
employ a non-perturbative approach by solving a non- 
linear equation for the electron Green function, which 
will allow us to ascertain the status of the previous re- 
sults obtained in perturbation theory and test our the- 
oretical predictions against several pieces of the existing 
experimental evidence. 

The semimetallic energy band structure of a single 
graphene sheet gives rise to the conduction and va- 
lence bands' touching each other in the two inequivalent 
ivf-points located at the corners of the hexagonal two- 
dimensional (2D) Brillouin zone. In the absence of in- 



teractions, the low-energy quasiparticle excitations with 
the momenta in the vicinity of these points labeled as 
i = 1,2 have linear dispersion Ep" 1 = ±vp, the velocity v 
being proportional to the width of the electronic 7r-band 
t re 2AeV f§. 

These excitations can be formally described by a pair of 
two-component (Weyl) spinors each carrying a spin 
index a, which are composed of the Bloch states residing 
on the two different sublattices of the bi-partite hexag- 
onal lattice of the graphene sheet. In what follows, we 
choose to combine them into one four-component Dirac 
spinor ty a = [^>\cr,ip2a) an d also treat the number of 
the spin components N as an adjustable parameter, the 
physical case corresponding to iV = 2. 

The use of the Dirac spinor representation allows 
one to cast the free quasiparticle Hamiltonian in the 
relativistic-like form where v is playing the role of the 
speed of light 
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where ^ a = ^^70 and the reducible representation of 
the 4x4 7-matrices 70,1,2 = (73, IT2, — iTi) <£> T3 given 
in terms of the triplet of the Pauli matrices Tj satis- 
fies the usual anticommutation relations: {TujTw} = 
2diag(l, -1,-1)1® 1. 

In the four-spinor representation, the electron 
Coulomb interaction reads as 
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Despite the apparent lack of the Lorentz invariance, 
Eq.(2) remains invariant under arbitrary U(2N) rota- 
tions of the 2iV-component vector La, ^Ra) composed 
of the chiral Dirac fermions defined as: ^l,r<x = |(1 ± 
75)^0-, where the matrix 75 = 1 ® T2 anticommutes with 
any %. 

The chiral invariance of Eqs.(l,2) brings about the pos- 
sibility of spontaneous chiral symmetry breaking (CSB), 
akin the phenomenon that has long been studied in 
the relativistic fermion theories. The CSB transition 
manifests itself in the appearance of a fermion mass 
and gapping of the fermion spectrum, thus breaking 
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the continuous chiral symmetry from U(2N) down to 
U{N) (g) U(N) and developing a non-zero expectation 

value < £^(r)vMr) >=< YZWMWM) ~ 
ip^ r (B)i(j (T (B)) >. The latter corresponds to the elec- 
tron density modulation which alternates between the 
two sublattices (A and B). 

In light of the above, one can identify the CSB order 
parameter with the site-centered CDW and thus relate it 
to the p = value of the gap function A p appearing in 
the renormalized (and, generally, non-Lorentz invariant) 
fermion Green function 
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where the interaction effects can also give rise to the non- 
trivial wave function (Z p ) and velocity (y p /v) renormal- 
ization factors. 

Because of its intrinsically non-perturbative nature, 
the phenomenon of CSB evades weak-coupling analysis 
based on perturbation theory. Nonetheless, similar to 
its relativistic counterpart ||, the CSB can be revealed 
by a non-perturbative solution of the system of non- 
linear equations for the fermion Green function (hereafter 
p = ej Q - vpj) 
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vertex function T p ^, and effective Coulomb interaction 
Vk = l/[(q/gv) + Nx(u> q.)] which gets strongly modified 
by the intra-layer polarization of the Dirac fermions 



Xk = Tr 
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A further analytical progress is hindered by the fact that, 
as a result of the interaction's Vk being explicitly non- 
Lorentz invariant, the gap function A p can feature sep- 
arate dependencies on the energy e and momentum p 
variables. Therefore, we choose to proceed directly with 
the finite temperature counterpart of Eq.(4), in which 
case the Lorentz invariance is broken regardless of the 
symmetry of the fermion interactions. 

In order to get a preliminary insight into the problem 
we resort to the same approximations as those made in 
the previous studies of CSB in the context of QED 3 [§. 
To this end, we first neglect the wave function, velocity, 
and vertex renormalizations (Z p — v p /v = r p .fc = 1) in 
Eq. (4) whose scalar part then becomes a closed equation 
for the fermion gap function. As shown in [pj , neglecting 
the above renormalizations in the gap equation suffices 
for establishing the existence of its non-trivial solution(s) 
and estimating a critical value N c of the only remaining 
free parameter, the number of fermion species. 

Taking the sum over the discrete Matsubara frequen- 
cies we then arrive at the momentum-dependent gap 
equation 
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where Ey, 



2 p 2 + A 2 . Next, we approximate the 

exact finite temperature fermion polarization (j^) by that 
computed in the massless case. By doing so, we overes- 
timate the contribution of the fermion momenta p < A p 
which is, however, unimportant, as long as the gap re- 
mains much smaller than the high-momentum cutoff A 
comparable to the maximum span of the Brillouin zone. 
As shown below, this condition is indeed satisfied for N 
close to the critical value N c . 

At A p = the fermion polarization is given by the 
approximate formula 

2T r 00 ^jq 
X(0,q) = — 2 dxln[2 cosh(— ^(l - x))] 
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which, for c = 161n2/7r, provides an up to a few percent 
accurate interpolation between the two opposite limits: 
vq 3> T where Eq.(Q) agrees with the zero-temperature 
result Xk °f V~k? and wq -C T where it exhibits thermal 
screening xo ccT || . 

Notably, at strong coupling (g ^> I) the screened 
Coulomb interaction Vk becomes independent of the bare 
coupling constant and assumes a universal form Vk ~ 
l/x(0,k) governed by the fermion polarization (7). 

Upon differentiating Eq.(6) with respect to the mo- 
mentum p one finds that for E p > T this non-linear 



(5) integral equation reduces to a linear differential one 
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which has to be supplemented by the boundary condi- 
tions A < oo and (A p + pdA p /dp)| p= A = 0. 

In turn, Eq.(||) can be readily identified with the radial 
Schroedinger equation for the s-wave zero energy level in 
the potential that behaves as oc 1/p 2 for p > T/v. From 
the textbook solution of this problem Q we infer that 
for N > N c = 8/ir there are two independent solutions 

A± cx i/pUi-sA--^/^)/ 2 , neither of which can satisfy 
the above boundary conditions. In contrast, for N < N c 
there exists a solution with infinitely many nodes, con- 
sistent with the infinite number of the negative energy 
levels in the I/p 2 -potential. Thus, only in this "center- 
ward downfall" regime does the solution 
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monotoneously decrease in the interval T/v < p < A and 
satisfy the boundary condition at p = A which reads as 
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where n is a positive integer. We note, in passing, that 
due to the formal similarity between the underlying equa- 
tions, a solution with the similar critical properties has 
also been discovered in the context of the 2D Cooper 
pairing near the antiferromagnetic instability [g . 

The highest critical temperature T C (N) below which 
the CSB order parameter sets in corresponds to n = 1: 
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According to Eq.(ll), at non-zero temperatures the crit- 
ical number of fermion species gets reduced: N c (0) — 
N c (T)^4ir 2 N c (0)/ln 2 (vA/T). 

In the case of QED3, the progressively more and more 
refined numerical simulations || have demonstrated that 
the results of the original analytical approach of Ref. 
H which yielded a solution similar to Eq.(9) remain 
robust against relaxing the above approximations and 
taking into account both the wave function renormaliza- 
tion and the vertex function satisfying the Ward identity 
r p ,o — Z p . 

Likewise, a numerical analysis of the coupled Eqs.(4) 
and (5) confirms the existence of the solution (9) in a 
whole domain bordered by the critical line (11) in the 
N-T plane @. 

Also, the numerically evaluated characteristic ratio 
2Ao/T c w 10 appears to be close to that found in Ref. || 
which is substantially greater than the BCS value corre- 
sponding to the solution A p = const of the gap equation 
with a momentum-independent kernel. 

As far as the nature of the CSB transition is concerned, 
the observed A^-dependence of the zero-momentum 
fermion gap Ao oc exp(—2ir/y /, N c /N — 1) prompts one to 
identify the breaking of the continuous chiral symmetry 
as a topological (Kosterlitz-Thouless-type) phase transi- 
tion in 2 + 1 dimensions [pj. Thus, the finite temperature 
CSB transition occurs between the two phases which are 
both chirally symmetrical, and therefore a bosonic Gold- 
stone mode must be present in the quasi-ordered phase. 

As regards the critical number of fermion species it- 
self, the recent symmetry-based argument made in the 
context of QED3 shows that the gap equation systemat- 
ically overestimates the actual value of N c which may, in 
fact, be as low as 3/2 [jl0| while the gap equation yields 

N QED = 32/7r 2 H 

The demonstrated formal relationship between the fi- 
nite temperature QED3 and the problem of a single 
graphene sheet suggests that in the latter case the ac- 
tual critical number of fermion species might also be less 
than two, hence no CSB occurs. 

Nonetheless, even at N > N c the nearby CSB transi- 
tion, albeit unreachable at any g, can still have a pro- 
found effect on the quasiparticle spectrum both above 



and below the crossover into the quantum-critical regime 
associated with the zero-temperature quantum-critical 
point at N c . 

In the quantum disordered (low-temperature) regime 
T < T*(N), where the crossover temperature T*(N) 
vanishes at A" — > A^+0 in the same manner as T C (N) 
given by Eq.(ll) for A" — * N c ~0, the only solution of 
Eqs.(4) and (5) is a massless fermion propagator which 
exhibits a suppression of the residue of the bare quasipar- 
ticle pole: Z p — > for p — > 0, while the velocity v p un- 
dergoes singular renormalization and monotonously in- 
creases with decreasing momentum, unlike in the Lorent- 
invariant QED3 where it remains constant. 

By contrast, in the quantum-critical regime T > 
T*(N) the fermion propagator features a simple pole, 
whereas the temperature- (but no longer momentum-) 
dependent factors Zt and vt control its residue and the 
effective velocity, respectively. The fermion damping is 
then determined by the self-consistent equation 
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which yields the universal solution: E(e) ~ max(e, T), in 
a general agreement with the time-resolved two-photon 
photoemission data taken in the energy range 0.4 < e < 
2eV [|J. One can expect that this characteristic signa- 
ture of the CSB-related quantum-critical behavior will be 
even more pronounced in the case of a graphite monolayer 
deposited on an insulating surface, whereas a conducting 
substrate would hamper the possibility of observing the 
linear damping due to strong metallic screening. 

The predicted linear damping should be possible to ob- 
serve in angular-resolved photoemission which can specif- 
ically probe the vicinity of the A'-points. On the con- 
trary, the angular-averaged data are going to be affected 
by such details of the graphite bandstructure as, e.g., 
the saddle point in the quasiparticle dispersion which 
occurs at e ~ 1.5eV if the momentum resides at one 
of the M-points of the Brillouin zone. This saddle 
point was recently argued to be a likely cause of the 
additional plateau-like feature observed in the angular- 
averaged E(e) jll| which is, therefore, unrelated to the 
many-body phenomena discussed in this Letter. 

In a stack of graphite layers with the inter-layer spacing 
d, the screened intra-layer Coulomb interaction remains 
dominated by the polarization \k only at q > 1/Ngd. 
At still lower momenta the kernel in Eq.(||) becomes less 
singular (Vfc oc 1/ ^/q) , thus reducing the range of tem- 
peratures and/or energies where the quantum-critical be- 
havior associated with the nearby CSB transition can be 
observed in electron photoemission. 

A finite inter-layer hopping t± ~ 0.27eV provides an- 
other cutoff below which the particle-hole pairing cor- 
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relations cease to drive the system towards the opening 
of the excitonic gap, and the quasiparticle damping be- 
comes quadratic in energy for e < t±. 

In contrast, the inter-layer Coulomb interaction has 
the opposite effect of nudging a stack of graphite layeres 
closer to the CSB instability. To elucidate this point, we 
recall that the common form of graphite has a crystal 
structure of well-separated hexagonal layers stacked in a 
staggered (ABAB . . .) configuration. As a result, each 
layer gets naturally divided into two sublattices formed 
by the atoms positioned just above and below the centers 
and corners of the hexagons in the two adjacent layers, 
respectively. 

Thus the inter-layer Coulomb repulsion strengthens 
the system's propensity towards developing the CDW in- 
stability by favoring spontaneous depletion of one of the 
two sublattices (accompanied by excess occupation of the 
complementary one) which alternates between the layers 
in order to keep the electrons in the neighboring layers 
as far apart as possible. 

Formally, one can incorporate this effect of the inter- 
layer Coulomb interaction into the effective single-layer 
description by adding a short-ranged four-fermion term 
HJ2a ^o-^ct) 2 into Eq.(2). With such a term present, 
the CSB transition occurs for any number of fermion 
species, including N > N c , provided that the strength of 
this coupling exceeds a certain critical value (A > X C (N)) 
which grows with N l| . 

By further elaborating on the solution of Eq.(6) ob- 
tained in the physical case N = 2 we also find out 
that, upon doping the system of graphite layers, the ex- 
citonic insulating ground state tends to spontaneously 
develop a non-zero spin polarization < xj^ a {A)\j} a {A) + 
^} 7 {B)^p l y{B) >~ fJ.S aao proportional to the chemical po- 
tential \i introduced by doping ||. 

This observation sheds light on the possible origin 
of weak ferromagnetism which was recently observed in 
highly oriented pyrolitic graphite (HOPG) [Q. Our 
fundings suggest that the latter might be not that dif- 
ferent from the mechanism proposed in the recent stud- 
ies of hexaborides believed to be 3D excitonic insulators 
p3[ . Notably, the authors of Ref. |T^] excluded mag- 
netic impurities as a possible cause of the ferromagnetic 
behavior of the magnetization hysteresis loops (also con- 
sistent with the electron spin resonance data) observed 
in the samples showing the insulator-like temperature de- 
pendence of the resistivity. Elaborating on the analysis of 
Rcf. jl3| we predict that if the excitonic instability proved 
to be at work, the weak ferromagnetism would have to 
disappear above a certain level of doping corresponding 
to the chemical potential /i c ~ Ao. 

In summary, we study the problem of the Coulomb 
interaction-driven electronic instabilities in layered 
graphite and propose a new explanation for the exper- 
imentally observed linear quasiparticle damping which 



might be a result of the relative proximity of a single 
graphene sheet to the zero-temperature quantum-critical 
point corresponding to the transition to the 2D exci- 
tonic insulator. In lightly doped layered graphite, the 
excitonic instability give rise to the formation of the 
site-centered CDW ground state exhibiting weak ferro- 
magnetism, as observed experimentally. Together with 
the recently proposed explanation [ fL4"[ of the apparent 
scmimetal-insulator transition in applied magnetic field 
]l5| as a phenomenon of the magnetic field-driven CSB, 
it lends a further support to the discovered formal re- 
lationship between the problem of layered graphite and 
the behaviors found in the relativistic theories of the 2D 
Dirac fermions. 
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